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Abstract 

Swimming of microorganisms is studied from a viewpoint of ex- 
tended objects (strings and membranes) swimming in the incompress- 
ible fluid of low Reynolds number. The flagellated motion is analyzed 
in two dimensional fluid, by using the method developed in the cili- 
ated motion with the Joukowski transformation. Discussion is given 
on the conserved charges and the algebra which are associated with the 
area (volume)- preserving diffeomorphisms giving the swimming mo- 
tion of microorganisms. It is also suggested that the iV-point string- 
and membrane-like amplitudes are useful for studying the collective 
swimming motion of microorganisms when fluctuation of the vorticity 
distribution exists in the sticky or rubber-like fluid. 
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1 Introduction 



In the high-energy physics, we usually study the world with the length scale 
of 10 -12 ~ 10 -29 /im, where field theories of particles as well as strings and 
sometimes exotic membranes play the important roles in describing theo- 
retically the phenomena of this tiny world. Unfortunately, we have never 
experimentally observed the dynamics of fission and fusion of strings and 
membranes, except for the old-aged hadronic strings. A lot of theoretical 
machineries have been so far developed to describe the dynamics of these 
extended objects, examples of which are conformal field theories based on Vi- 
rasoro, Kac-Moody and W-algebras (area-preserving diffeomorphisms) , the 
field theory of strings and so on. Therefore, it is by no means a waste of time 
to search for another region of physics where the knowledge of high-energy 
physics will probably be useful. 

In this paper, we take up the swimming problem of microorganisms as 
such an interesting application, where we know the following fact on the 
ways of swimming |l] : There exist only three different universality classes of 
the swimming ways of microorganisms; (1) Swimming with cilia is adopted 
by the spherical organisms with the length scale of 20 ~ 2 x 10 4 /mi, an 
example of which is Paramecium; (2) the smaller microorganisms with the 
size of 1 ~ 50/zm swim with flagella, an example of which is the sperm; 
(3) the bacteria with the size of 0.2 ~ 5/mi swim with bacterial flagella 
the motion of which resembles the screwing of the wine-opener. Why is it 
possible that such simple classification be realized in the swimming problem 
of microorganisms? 

Recently, Shapere and Wilczek have studied the swimming of mi- 
croorganisms from the gauge theoretical view point; Knowing one cycle of 
the swimming motion gives a closed path on the space of shapes, they have 
found that the net self-propulsion of translation and rotation in the low 
Reynolds liquid are expressed in terms of the Wilson operator on the closed 
path. Starting from their works, we will develop the ciliated and flagel- 
lated motions and study on our problem from the view point of algebraic 
structure existing in the deformation operation of the microorganisms in the 
low Reynolds liquid, as well as from the swimming dynamics of a group of 
mocroorganisms, which may be connected to the iV-point correlation func- 
tion of string and membrane theories. 

Let us recall the hydrodynamics. The Reynolds number R, giving the 
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ratio of the kinetic term over the viscosity term by R = pVL/p, where 
p, V, L, and p denote density, typical velocity, typical length scale and the 
coefficient of viscosity, respectively, satisfies 

R < 1 x [L/mm], (1) 

for the velocity of microorganisms V < lmm/s with p = lg/cm 3 and p = 
lCT 3 Pas of water. Therefore, for the microorganisms with the length L <C 
1mm, we have R <C 1 (the low Reynolds number fluid) which leads to the 
following equations of motion for the incompressible fluid: 

V • v = 0, (2) 

and 

Av = -Vp, (3) 
p 

or equivalently 

A(V x v) = 0, (4) 

where p is the pressure and v(x) is the velocity field of the fluid. The surface 
of a microorganism swimming in the D = 3 dimensional fluid of the real 
world forms a closed membrane at a fixed times t, the position of which can 
be parametrized by introducing (D — 1) parameters £ l (i = 1, • • • , D — 1) as 
X M = X^it; • • • , ^' D_1 ). It is sometimes instructive to consider the D = 2 
dimensional fluid. Then, the surface of a ciliate (flagellate) becomes a closed 
(open) string and its position can be described by a complex number 

Z = x 1 +ix 2 = Z(t;9), 

with —n < 9 < it. In the sticky fluid of R = 0, there is no slipping between 
the surface of a microorganisim and the fluid, namely, we have the matching 
condition 

v(x = X(f;0) = X(f;0, (5) 
or in the general coordinate system with the metric tensor g^ v (x) 

g^(x)v,(x)\ x=x(W =X»(t;0- (6) 

Now we will exemplify the ciliated as well as flagellated motion of mi- 
croorganisms in D = 2 dimensional fluid. The ciliated motion was attacked 
by Shapere and Wilczek several years ago, so that our main concern will 
be the flagellate motion, but we first discuss the ciliated motion from a bit 
different viewpoint by summarizing our notations. 
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2 The ciliated motion in the D=2 fluid 



Time evolution of the envelope of the cilia can be viewed as a small but 
time-dependent deformation of a unit circle in a properly chosen scale, 

Z(t,9) = s + a(t,s), (7) 

where s = e %e and a(t, s) is arbitrary temporally periodic function with period 
T satisfying \a(t, s)\ <C 1 with — n < 9 < ir. The complex representation of 
the covariant velocity vector u M can be denoted as 

2v E (z,z) = (vx + w 2 ){z,z) (8) 
2v z (z,z) = (vi-iv 2 ){z,z), (9) 

since their transformation property under the conformal transformation z = 
f(w) reads 

Ow 

v s (z,z) = —Vyjiw.w) (10) 

v z (z,z) = ^-v w (w,w). (11) 
oz 

The matching condition @ is then given by 

2v g (s + a(t, s), s' 1 + a(t, sfj = a(t, s). (12) 



The factors 2 in Eqs.(^|), (|9p and ( Jl^ ) come from the flat two dimensional 
metric expressed in terms of z and z, that is, g zz (z,z) = 2. Perturbatively, 
Eq.(fl2|) can be solved as 

2v- z {s) = a(t, s) - 2 \a(t, s)d z v^ + ^[€s)d- z v z {0) ] (13) 

J z=s 

with the solution of the lowest order v z ^(z,z) satisfying 

2v^°\z,z)\ z=s = a(t,s). (14) 

The general solution of the velocity vector for the R ~ 0, incompressible 
liquid is already known; the equations of motion @ and (fl) are then 

d z v z + d z v z = (15) 
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and 

kd z d z v- z = -d^ (16) 
/i 

or equivalently 

d z d- z (d z v- z -d 2 v z ) = 0. (17) 

Introducing the stream function U(z,z), Eq.fll5|) is solved as v z = d z U and 
v z = —d z U, and Eq. (pj) in terms of [/, (d z d z ) 2 U(z, z) = 0, has a general 
solution with two arbitrary holomorphic functions <f)\(z) and <fro(z) as 



U(z, z) = zfa(z) ~ z<f>i{z) + <f>o(z) - <f>o(z). (18) 
Now, we have the general solution for v z as 



v s (z, z) = faiz) - z4J x {z) + fa(z), (19) 

where <j>2{z) = <fi' (z). If we require the velocity field v z be finite at spacial 
infinity, we have the expansion, 

0i(z) = £ a fc ^ +1 and 2 (z) = £ 6^ fe+1 . (20) 

fc<0 fc<-l 

Following Shapere and Wilczek, self-propulsion of translation and rotation 
against the fluid can be extracted as the counterflows of — a_i and — Im6_2, 
since the translational and rotational flows, a_i and Im6_ 2 , surviving even 
at spacial infinity after solving Eq.([5|) or (12) should be cancelled by the self- 
propulsion of the microorganism. Here we will find the following useful prop- 
erty; if the velocity vector is known on the unit circle as v(s) = Y^t=-oo v kS k+1 , 
then the velocity vector at any place can be easily determined as 



v g (z, z) = v { ~\z) + v i+ \z) + (z- 1 - z)i/(-)(z), (21) 
with the help of 



v y '{z 



Y.VkZ'*- 1 (22) 



fc>0 



v (-)( z ) = J2v k z k+1 . (23) 

fc<0 
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Then from Eqs.([HJ), (0), and (|21|), we obtain a general expression of v s 
caused by the time dependent small deformation of a(t, s) as 

2v 2 (s) = a(t,s) 

- {a{t, s) {d a 6^-\t, s) - d a a(-)(t,s)) 

+ ajt^s) (d- s a (+) (t, s) - s 2 d s a^(t,s)) } . (24) 

From this expression, the net translationally swimming velocity v^ lha ^ of the 
ciliated microorganism reads 

2 4 cilia) = -a (t) 

^l(a n a — n+1 C^n^n—l Ot n Ol.— n-\s) 

- nd n o^~[, (25) 

n>l 

where a n {t) is defined by a(t, s) = Y.n=-co a nit)s n . On the other hand, the 
net angular momentum t>^ lba ^ gained by the microorganism from the fluid 
becomes 

Im {ai(t) 

- Y n(a n a^ n+2 - d n a n - a n a>- n+2 ) 

n<l 

+ nci nC^ \ • (26) 
n>l J 

The net translation and rotation resulted after the period T come from 0(a 2 ) 
terms since the 0(a) terms cancell after the time integration over the period. 

3 The flagellated motion in D=2 fluid 

Microorganisms swimming using a single flagellum can be viewed as an open 
string with two end points, H and T, where H and T represent the head 
and the tail-end of a flagellum, respectively. Our discussion will be given 
by assuming that the distance between H and T is time-independent and is 
chosen to be 4 in a proper length scale. This assumption can be shown to be 



2v 



(cilia) 
R 
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valid for the flagellated motion by small deformations in the incompressible 
fluid 0. Then, at anytime t, we can take a complex plane of z, where H and 
T are fixed on z = 2 and —2, respectively. This coordinate system z can be 
viewed as that of the space of standard shapes of Shapere and Wilczek. Time 
dependent, but small deformation of the flagellate can be parametrized as 

Z(t,0) = 2(cos9 + isin9a(t,9)), (27) 

where the small deformation a(t, 9) can be taken to be a real number[| satis- 
fying 

a(t,9) = -a(t,-9). (28) 

Here, we parametrize the position of the flagellum twice, starting from the 
end point T at 9 = —n, coming to the head H at 9 = 0, and returning to T 
again at 9 = n. Motion of the two branches corresponding to — n < 9 < and 



7r > 9 > should move coincidentally, which requires the condition (|28f) . The 
Joukowski transformation z = z(w) = w + w -1 , separates the two coincident 
branches in the z plane to form lower and upper parts of an unit circle in 
the w plane, outside domain of which we are able to study the swimming 
problem of the flagellate in a quite similar fashion to that of the ciliate, 
The parametrization of our microorganism in the w plane corresponding to 



Eq . (|27|) is now 

W{t,9) =e ie (l + a(t,9)) + 0(a 2 ). (29) 

The general solution of the velocity field Vg in the z plane can be transferred 
to that of in the w plane, by using the transformation property fllCf ), that 
is, we have generally 

dw 

Vw{w,w) = i)i{w) - z(w)—ijj[(w) + if) 2 {w), (30) 

where the two arbitrary holomorphic functions have the following expansion 
Mw) = E a ^ k+1 and ^ 2 {w) = £ b k w k+1 . (31) 

fc<0 fc<-l 



1 When a is taken to be a complex number, the length of the flagellum is locally change- 
able at 0(a). For such an elastic flagellum, we have similar results to that of the ciliated 
motion. In case of real a, its length is locally preserved at 0(a), giving a non-elastic 
flagellum, which is the more realistic one. 
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In the study of the ciliated motion, the property ( |2T| ) was extremely use- 
ful; the corresponding expression in the flagellated motion can also be ob- 
tained from a given velocity vector on the unit circle in the w plane = 
J2t=-oo v kS k+1 , by using the general solution of (|30|) and (|31|). In terms of 
velocity field of the real space, we have 

Vz{w,w) (32) 
= (1 - w 2 )- 1 (a + pw- 1 - w~ 2 v^\w)) 

+ (1 - w- 2 )- 1 (aw- 2 + (3m- 1 + v i+) (w)) 
+ (w - w _1 ) _3 {(w + w' 1 ) -(w + w' 1 )} 

x [(3(w + w' 1 ) + 2a- 2v(-\w) + (w - ttT 1 )^-)'^)} , (33) 
where v'^ are defined by 



(+)(^ = J2v k w~ k - 1 (34) 



k>0 



v y '[w 



Y,VkW k+1 . (35) 



k<0 



The constants a and (3 in Eq. (|33|) are fixed so that v 2 has no singularities on 
the special points of the head H and tail-end T. From this requirement, we 
have 



1 

a = - I v 



,(-)(!) +„(-)(_!)) (36) 

KO 

P = \^-\l)-v^-\-l))=Y J V2u (37) 

KO 



as well as 



Vz(w,w)\ w= u,=±i = 0, (38) 



which is quite consistent to our requirement of fixing the head and tail-end 
position in z and w planes. 

Next task is to solve the matching condition of @ which is written in the 
z plane as 

2v- z (W(t,6),W(t,9)) = Z(t,9). (39) 
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Perturbatively, the solution of p5|) reads 
2v- z {s) = (s-s _1 )d(t,0) 



sa(t, 9)d w vf^ + s 9)dwV 



,(°) 



with the help of the lowest order solution f^^w, w) satisfying 



2v 



(0)| 



(40) 



(41) 



on the unit circle (s = e td ). Using the mode expansion satisfying Eq. 



a 



(t, 9) — a n (t) sin n6*, 



n=l 



the lowest order solution is obtained; 



2iVg (w, w) 



(42) 



(43) 



n>l n>l 

-1\ ^T. I „-T.-l> 



+{(w + iy x ) - (w + u) L )} 
w + w" 1 



E 

n>l 



x < > na„iu - 



w — w 



(44) 



n>l 



Then, from Eqs.([40"|)^([44|), we have derived the velocity of the fluid which 
incorporates perturbatively the small but time- dependent deformation of the 
flagellate as 



2v- 7 



2i sin 9 V] a n (t) sin n9 



n>l 



—4 ^ a n (t)d m (£) sinn#(cos#sinm# + msin# cosm#). (45) 

n,m,>l 

Now, we are able to determine the net swimming velocity y T ( fla s ella ) gained 
by the flagellate motion of microorganisms: 



2 t ,C[ ia s ella ) 



-iati - ma m a m+ i + ^ ma m « m _i, 

m>l m>2 



(46) 
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On the other hand, the angular momentum u R C aa s ella ) is given by 

Wflagella) = ~OL 2 . (47) 

After the time integration over the period T, y R ( fla s ella ) vanishes since in 
our first order approximation, the length of the flagellum is fixed in the 
incompressible fluid. Therefore the second order approximation is necessary 
for the non- vanishing y R ( fla s ella ). 



4 The selection rules and the symmetry of 
microorganisms' swimming 



Even though the results in Eqns.(P5|), (|2q), (fiq) and (fTfl ) are obtained per- 
turbatively, we are able to read from them the characteristics of the microor- 
ganisms' swimming; In order for the ciliates to swim or rotate, they need the 
coexistence of the two different Fourier modes of ni, n 2 . The selection rules 
for the allowed (n 1 ,n 2 ) conbinations are 

f) \n\ — n 2 \ — 1 or |m + n 2 — 2\ = 1 for the cililated translating) 
if) n\ = n 2 or \n\ + n 2 \ = 2 for the ciliated rotation, (49) 

The corresponding selection rules for the flagellate motion are 

in) |rii — 77.2] = 1 for the flagellate translation, (50) 

where the Fourier modes are sinn# in this case. 

Viewing these selection rules, we are tempted to elucidate the algebraic 
structure possibly existing in the background of the swimming mechanism. It 
is similar to the Virasoro algebra, but is different from it. For such a purpose, 
introduction of the "action" will be convenient. The "action" S reproducing 
the classical equations of motion of the swimming of iV microorganisms in 
the incompressible liquid with low Reynolds number may be given by 

N 

s N = ^y^y^-^p^^^ 



i=l 



1 

Ji y "' ' 4 



p(x)d u v u (x) + -u^(x)uu^(x) 



(51) 
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where the velocity field oj^x) is given by 

W/tvix) = d^v v - (52) 

whose D = 2 expression is uj z ,z{z, z) = d z v z — d z v z . 

We have introduced the parameter a' so as to make S dimensionless, 
where a' has the dimension of (Length) D (Time) -2 . The reason why we have 
used the notation a', familiar in the string theories to describe the Regge 
slope, will be understood later. The Lagrange multiplier fields of Pfy(t;£(i)) 
(i = 1,---,N, /i = 1, ■•■,D) guarantee the matching condition of or 
@ for i-th microorganism at any time, and the pressure p(x) is also such 
multiplier giving the incompressibility given in Eq.(|2|). The field equation 
(H) can be easily reproduced. In the action ([51]), time t appears only in the 
first term of representing the matching conditions, that is, the time evolution 
is triggerd only by the self-motion of the microorganisms, of which influence 
spreads instantaneously over the whole space and causes the change of the 
fluid velocity there. Because of the lacking of the kinetic term, we may 



call S as the "action". The additional metric contribution such as y g{x) is 
only relevant for the curved space, an example of which has appeared in the 
flagellate swimming on the w plane. The later discussion is given for the flat 
metric. 

Now, we will define the following local transformation at a fixed time t: 

SXfo = A"(X w (f; £(,■))) (53) 

5P® = (54) 

5v»(x) = \»(x) (55) 

6p(x) = k(x), (56) 

where we have assumed that the transformation parameters A M (x) and k(x) 
are restricted by the equations of motions, 

8^{x) = 0, and dpWvix) = -8 v k(x), (57) 

where = d^Xy — d u X^ is the vorticity for A M . Meaning of the transfor- 



mations (|53|) ~ (^6|) are quite simple; the deformation of the microorganisms 
d53|) triggers the increase of the velocity field fl55|) and of pressure (|56|) so 
that they can be consistent with the incompressible fluid dynamics of the 
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low Reynolds number. It is also important to note that the succession of 
these time-independent transformations result in the time evolution of our 
problem. Therefore, the transformations (^)~(^) resemble the ordinary 
canonical transformation generated by the Hamiltonian. 

Now, let us find the conserved current associated with the transformations 



D~dHD- ^ we denote the Lagrangian of the fluid itself as £/; 

C f = — pd^ + -u^ v , 

fi 4 



then we have for the transformations 



5C f = 



(58) 



(59) 



Following the usual method, the conserved current associated with the pa- 
rameters X tl (x) and k(x) is obtained: 



5C t 



-5v h 



A 
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1 



<Tp )K - A 



(60) 



We should notice that the conservation of this current can be directly proved 
with the help of the equations of motion (Eqs. (0) and (J3|)), namely, 



d^v„ = 



(61) 



and 



0. 



(62) 



as well as Eq. (p7|) . Instead of Eq.(J62[), we can adopt the conservation of the 
stress tensor 

(63) 



d»T>" = 0, 



with 



= fi(dpV v + d v v^) - g^p. (64) 

Therefore, the conservation of the current J m [ k ,a"](^) is essentially the reflec- 
tion of the incompressibility and the conservation of the stress tensor or the 
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balancing of the stress force of the fluid. Using the current conservation in the 
outer region of the microorganisms, Gauss theorem leads to a conservation 
law; 

N 

£%) M " ] = 0, (65) 

i=l 

where 

Q {i) [K ' X " ] =l dP-^J^ix) (66) 

is the integration over the surface of the z-th microorganism. It should 
be cautious about the flow at spacial infinity which corresponds to the net 
translation and rotation of the group motion of the microorganisms. To 
take into account of this effect, it is convenient to include an imaginary 
microorganism, say % = N, at spacial infinity: Stm may be an infinity large 
envelope of the whole space, or in the compactified space it can be an ordinary 
microorganism located at the infinity point. 

The transformation also generates the deformation of the shapes of mi- 
croorganism fl53p. The generator of this deformation per unit time can be 
written as 

V] = /^WfeO)j^ir (67) 

which gives the volume (area for D = 2) preserving diffeomorphisms ow- 



ing solely to the incompressibility condition in Eq.(|57|). The second con- 
dition in Eq.flSTD adds the further restriction on L\\vy. By the help of the 
stream function <J\(x), the incompressibility condition is automatically sat- 
isfied through A M = e^ uX d u cr\, so that the Eq.(^) becomes the constraint on 
the stream function 

A(Ag^ - d lt d v )a^{x) = 0, (68) 
or in D = 2, in terms of the only non-vanishing component a = o~3 

A 2 a(x) = 0. (69) 

This constraint has been already solved generally in D = 3 and D = 2 fluid. 

Here we will discuss a little more on the nature of the charge Q^ K ' X "^ 
in Eq.(|^) as well as the algebraic structure of the generator in Eq . (|67|) , 
taking up the simple case of D=2. Then, we can use the knowledge given in 
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Eqs.(|H)|)^(|i~9|). The current J M (x) given in Eq.(|60|) reads in D=2 



J"[k,\](z, z) = 2[ 2u Sz p ) X z - 2 2X 2z « 



(70) 



where the factors of 2 come from the metric g zz = 2 needed for raising and 
lowering the indices. The general expression of the parameters k and A are 
obtained as 



X z = e x (z) -ze 1 '(z) + e 2 (z), 



Ai 



2{e 1 '{z)-e l \z)) 



and 



k = -Afi{e 1 '{z) + e 1 '(z)), 
where we take the following most general expansion 



+oo 



ei(z) = d L \nz + ^ c k z 

k=—oo 



k+l 



and 



+oo 



(71) 
(72) 

(73) 
(74) 
(75) 



e 2 (z) = d L \nz + d kZ k+1 , 

k=—oo 

admitting the locally-defined but single-valued transformations. The fields 
v z , co>zz, and p have of course the expressions similar to (|?l~l)^(|7"3"l)'s, where e± 
and e 2 are replaced by 0i and (p 2 in Eq. (PU[) . Now we have 



■J Z [eue 2 ](z,z) = 16 



(76) 



The corresponding charge defined in Eq.(^) on the circle, \z\ = r, takes the 
following form; 



n (\z\=r) 



dz rz r dz , 
J + 4> — J 



=r 2ni J\ z \=r 2m 
= %{-(a^id L + a~[d L ) 

+ (k + l)(a fc rf_ fc _ 2 + + hc-k-2 + hcZ^Ji)}. (77) 



fc<-i 
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This means that there are infinite number of conserved charges which are just 
equal to the coefficients {cik}k<o and {bk}k<-i of </>i(z) and <f>2{z). In other 
words, two functions (f>i(z) and <fi 2 (z) give a solution of the fluid dynamics in 
the D = 2 domain, so that if r is considered as the evolution parameter, the 
functional form of them is conserved. 

Next, we will study the algebraic structure of the operator in Eq. (|67|) . 
Let's introduce the local version of its operator 

L M ee A"(aO^, (78) 

which becomes in D = 2 as 

L a = 2(X 2 d z + \ z d z ) 

= 2(d- z ad z -d z ad- z ) (79) 

when it is written in terms of the stream function a(z, z) having the follwing 
general form: 



rz rz 

a(z, z) = Z€\ — ze\{z) + / dze 2 (z) — / dze 2 (z). (80) 

Then, L a can be considered as a Liouville operator of a dynamical system 
moving in the phase space of (z, z), having —a(z, z) as its Hamiltonian. In- 
variance of the phase volume during the temporal evolution of the dynamical 
system shows that L a is really the area preserving diffeomorphism. The com- 
mutation relation is simple, namely 

[L ai , L a , 2 ] = —L{ C1)(J2 }, (81) 

where the {a Xl cr 2 } is the Poisson bracket defined by 

{fi, 0-2} = d z axd z a 2 - d z a x d z o 2 - (82) 

The characteristic of our area-preserving diffeomorphism is summarized in 
Eq.flgQD which can be expanded as 

<t(z,z) = Y, c ^z k+l -J2^zz k+1 

k k 

+d- 2 lnz — d_ 2 \nz + d^zihizz — 1) — d^zi^uzz — 1) 

+ E 7^^ +2 - E j^2 dkzk+2 - (83) 
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The generators corresponding to each term of the expansion may be denoted 

as 

L Ck = L z k+i s , Lj-^=L\ n2 , and L^ L — L z inz2-i, etc (84) 

Most of the generators are associated with the monomial z mi z m2 which can 
be denoted as Lr mi)Tn2 y Then, the following commutation relation is easily 
understood. 

[i(mi,m 2 )) ^(ni,n 2 )] = — ( m l n 2 ~ m 2 n l) L( mi+ril -i i m2+ „ 2 _i) , (85) 

which defines the so-called VFi+oo-algebra, having L(o, o) as its U(l) operator. 
The overall translations and rotation are generated by L C l , L 5 _ 1 and Lfe_ 2 + 
Li_ 2 , respectively. Virasoro-like algebras are included among the generators; 

k. L cJ = -(m-n)L Cm+n , (86) 

l L c m , L 5n ] = (m-n)L 5m+n , (87) 

but 

[L Cm , L 5n ] = -(mn + m + n)L (m+ljn+1) (88) 

results in the appearance of a new generator L( m+ i i n+ i) . It is quite inter- 
esting to understand that net translation and rotation for the ciliated mi- 
croorganisms are generated by the generator £( m +i, n +i), where |n ±m| = 1 
corresponds to translation ( |49D and n = m or \n + m| = 2 corresponds 
to rotation [0]. As yet, the role of the generators {L a } in the swimming 
problem is not manifestly clear but the algebra can depend on the shape of 
the microorganisms that really happens in the case of area-preserving diffeo- 
morphisms a part of which is included in our algebra. We could guess that 
the study on the classification of possible area (volume)-preserving algebras 
of this kind may lead to the understanding of three universality classes of 
microorganisms' swimming discussed in the beginning. 

5 Collective motion of microorganisms and 
N-point correlation function 

Finally, starting from the "action" ([51]) of a group of microorganisms, we will 
study a possibility of using iV-point correlation functions similar to those 
of string and menbrane in the study of swimming motion in the D = 2 
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and D = 3 fluid, respectively. Consider the situation where the vortices 
are created and annihilated, so that the probability of having the vortex 
distribution (field) uj^ v {x) is given by 

P[^(x)} ~ exp {-^7 / d D x^(xj iuv^K^)} , (89) 

where ia! is the external parameter controlling the fluctuation of the vortex 
distribution. [The a 1 — > limit corresponds to the classical limit without the 
fluctuation.] In this situation we should sum over all the possible configura- 
tions of the velocity fields with Eq. (^3) as their probability. The probability of 
having N microorganisms with their surfaces located at X^,X^ 2 ), ■ ■ ■ ,X( N y 
and with their time derivatives -XVi), Xm, ■ ■ ■ , Xrm, is given by the following 
iV-point correlation function (amplitude) 

Gn(X(i), X(i); X(2),X( 2 ); X^,X^n)) = Z n /Z . (90) 

Here the partition function Z N is defined by 

Z N = J VP®(t,Z {i) ) J Vp(x) J Vv^x) exp{iS N }, (91) 

with the "action" Sn in Eq-(H)- In the path integrations over P^(t; £(j)) and 
p(x) in Eq.(|9T|) guarantee the matching condition (|6]) and the incompress- 
ibility of the fluid, respectively. Moving to the momentum representation 
natually, we have 

Gn(X(i), X^y, X( N y X^) 

= J VPf exp j*f] Jdt J dP-^P^X^ 
xG N (X (1) ,P^; ... ;X {N) ,PW), (92) 

where 

G N (X {lh pW;--.;X {N) ,pW) 
= Z l f Vv"{x)\{8{dv{x)) 
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If 1 

x exp / d D x -v^Av u (93) 

P 27ria' J 2, 

= exp 



2 W x ~J2 J dt (i) J d D ~^ (t) J dt {j) J d ' 1 ^ 



x P«(*( ; - ^(to; M • ( 94 ) 

Here the Green'n function G^{x — x') satisfies 

AG^(i — x') = 8^ D \x — x'), (95) 

and the transverse condition reflecting the incompressibility of the fluid, 
namely 

d^(x - x') = & v G^{x - x') = 0. (96) 

Now, we have obtained iV-point correlation function for the collective swim- 
ming of N microorganisms. It is quite similar to the iV-point function of 
strings || for D = 2 case and membranes || for D = 3 case. It is also 
related to the string field theory, since the incoming and outcoming strings 
are not point-like, but the Reggeons. More precisely, refinement is necessary 
on the treatment of functional measure and boundary condition in Eq. (|9"3"|) . 
There is a possibility that the boundary condition terms changes Ag^ into 
Ag^ v + Svdy, - g^s x d x with 

N 

i=i J 

It is also interesting to consider the multi-loop amplitudes if the meaning of 
handle in the space of velocity field could be well understood. 

Here we will discuss the use of iV-point correlation function (j9^) in the 
collective swimming motion of microorganisms. If G^ represents the prob- 
ability of having N microorganisms whose surfaces are located at Jm with 
velocity (i — 1, • - • , N), then it can be viewed as the probability distri- 
bution of Xtv of the imaginary microorganism N located at spacial infinity 
under the given data of X { and X^i — 1, ■ ■ ■ , N — 1). Following the usual 
strategy, the counterflow — Xn can be indentified to the collective swimming 
motion of N — 1 microorganisms. Therefore, the averaged collective swim- 
ming motion over the fluctuation distribution is given by 

— (Xn) — ~ X/ -^nGn{Xi, ' " ' ; Xn, X n ). (97) 
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It is also interesting to study the response of the Appoint correlation function 
Gn, the Ward indentity, for the deformation of the shapes of microorganisms 
in Eq.O—O's. 

We hope much progress will come out along the strategies proposed in 
this article. 
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